The fate of bound systems in phantom and quintessence cosmologies 
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We study analytically and numerically the evolution of bound systems in universes with accel- 
erating expansion where the acceleration either increases with time towards a Big Rip singularity 
(phantom cosmologies) or decreases with time (quintessence). We confirm the finding of Caldwell 
et. al. U that bound structures get dissociated in phantom cosmologies but we demonstrate that 
this happens earlier than anticipated in Ref. Q. In particular we find that the 'rip time' when 
a bound system gets unbounded is not the time when the repulsive phantom energy gravitational 
potential due to the average (p + 3p) balances the attractive gravitational potential of the mass M 
of the system. Instead, the 'rip time' is the time when the minimum of the time dependent effec- 
tive potential (including the centrifugal term) disappears. For the Milky Way galaxy this happens 
approximately 180Myrs before the Big Rip singularity instead of approximately 60Myrs indicated 
in [J for a phantom cosmology with w=-1.5. A numerical reconstruction of the dissociating bound 
orbits is presented. 

PACS numbers: 



I. INTRODUCTION 

There is mounting observational evidence that the uni- 
verse has entered a phase of accelerating expansion (the 
scale factor obeys a > 0) and that the total amount of 
clustered matter in the universe is not sufficient for its 
small average spatial curvature. This converging obser- 
vational evidence comes from a diverse set of cosmolog- 
ical data which includes observations of type la super- 
novae [2j, large scale redshift surveys [3j and measure- 
ments of the cosmic microwave background (CMB) tem- 
perature fluctuations spectrum Q. The observed accel- 
erating expansion and flatness of the universe, requires 
either a modified theory of gravity^ or, in the context 
of standard general relativity, the existence of a smooth 
energy component with negative pressure termed 'dark 
energy '[(|. This component is usually described by an 
equation of state parameter w = ^ (the ratio of the ho- 
mogeneous dark energy pressure p over the energy den- 
sity p). For cosmic acceleration, a value of w < — | is 
required as indicated by the Friedmann equation 



4ttG 



(P + 3p) 



(1.1) 



Current observational bounds 0, on the value of the 
dark energy equation of state parameter itf(to) at the 
present time to yield 



-1.48 < w(t ) < -0.72 



(1.2) 



at the 95% confidence level. The role of dark energy 
can be played by any physical field with positive energy 
and negative pressure which violates the strong energy 
condition p + 3p > (w > — |). 
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Quintessence scalar fields^] (—1 < w < — |) violate 
the strong energy condition but not the dominant en- 
ergy condition p + p > 0. Their energy density scales 
down with the cosmic expansion and so does the cosmic 
acceleration rate. Phantom fields 0,^1^1 (w < — 1) vi- 
olate the strong energy condition, the dominated energy 
condition and maybe physically unstable. However, they 
are also consistent with current cosmological data and 
according to recent studiesQ they are favored over their 
quintessence counterparts. In contrast to quintessence 
fields, the energy density of phantom fields increases with 
time and so does the predicted expansion acceleration 
rate -. This monotonically increasing acceleration rate 
of the expansion may be shown to lead to a novel kind 
of singularity which occurs at a finite future time and 
is characterized by divergences of the scale factor a, the 
Hubble parameter H its derivative H and the scalar cur- 
vature. This singularity has been called 'Big Smash' ^| 
the first time it was discussed and 'Big Rip' (here- 
after CKW) in a more recent study. Even though there 
are mechanisms by which the 'Big Rip' singularity could 
be avoided pH] it remains an interesting possible fate of 
the universe. 

An immediate consequence of the very rapid expan- 
sion rate as the Big Rip singularity is approached is the 
dissociation of bound systems due to the buildup of re- 
pulsive negative pressure in the interior of these systems. 
This observation was first made in where a qualitative 
study of the dissociation times for various bound systems 
was also made. 

The quantitative study however of the evolution of a 
bound system in an expanding universe remains an issue 
of current research. In particular the question of whether 
the expansion of the universe affects gravitationally 
bound systems like clusters, galaxies or planetary sys- 
tcms has been addressed in several early (fl El Q |l7| 
and recent 0, 0, |2(j studies. The recurrent attention 
paid to this issue indicates that to this point a definitive 
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answer is still lacking. The prevalent perception however 
is that the physics of systems which are small compared 
to the radius of curvature of the cosmological background 
is essentially unaffected by the expansion of the universe. 
Even though this perception is applicable in the case of 
phantom cosmology, it is not applicable as the Big Rip 
singularity is approached and curvature increases rapidly. 

A qualitative approach to the problem of bound sys- 
tem dissociation in phantom cosmologies was made by 
CKW where it was assumed that a bound system be- 
comes unbound when the source of the repulsive phan- 
tom energy gravitational potential for a bound system of 
size R (— |7r(p + 3p)R 3 ) balances the attractive source of 
the ordinary gravitational potential (the mass M of the 
system). Even though this assumption is qualitatively 
useful as it provides a physical understanding as to why 
does a bound system dissociate due to phantom energy 
repulsion, quantitatively it leads to incorrect estimates 
of the dissociation times because it ignores the effects of 
the centrifugal barrier in the effective potential of bound 
systems. A quantitative analysis should make use of the 
particle equations of motion in the local inertial frame 
based on the geodesic deviation equation. Using such an 
equation, the time dependent effective potential of two 
body bound systems may be shown to consist of three 
terms, the static attractive gravitational mass term, the 
static repulsive centrifugal term and the time-dependent 
repulsive dark energy term. The stability of a bound sys- 
tem depends on the existence of a minimum for the above 
effective potential. In the case of quintessence (w > —1) 
the minimum of the effective potential can not disappear 
at any future time because the repulsive time-dependent 
term decreases with time. For phantom energy however 
(w < —1) the repulsive time-dependent term increases 
with time and at some critical time dominates over the 
other two terms and destroys the minimum of the effec- 
tive potential dissociating at the same time the corre- 
sponding bound system. It should be stressed that this 
critical time when the effective potential minimum dis- 
appears is not the time when the repulsive dark energy 
term balances the attractive gravitational mass term as 
assumed in CKW. This will be demonstrated in detail, 
analytically and numerically in the following sections. 

The structure of this paper is the following: In section 
2 we use a metric that interpolates between the static 
Schwarzschild at small scales and a general time depen- 
dent Friedmann universe metric at large scales to derive 
the geodesies of a test particle in the Newtonian limit. 
We then focus on the particular Friedmann universe that 
contains a two component cosmic fluid (matter and dark 
energy) and derive the form of the scale factor for any 
value of constant w (the equation of state parameter). 
Using this form of the scale factor in the previously de- 
rived geodesic equation we find the equation of motion for 
a two body bound system in an expanding quintessence 
or phantom cosmological setup. In section 3 we study 
this equation of motion and derive the evolution of the 
radius of two body bound systems with distance and mass 



scales corresponding to the Solar System, the Milky Way 
Galaxy and the Coma Cluster. We also derive analyti- 
cally the dissociation time as a function of w for a bound 
system of given mass and radius and test our result using 
numerical simulations of the above three types of bound 
systems (Solar System, Milky Way and Coma Cluster). 
Our analytical result for the time difference 'Big Rip time 
- Dissociation time' differs from the corresponding result 
of CKW by a factor of 3.08 (we find the dissociation 
to occur earlier). The time evolution of the bound sys- 
tem effective potential is also shown. Our results are 
extended to the case of planar mass sources and to the 
case of quintessence (w > — 1). In the later case it is 
shown that the size change of the bound system due to 
the expansion is negligible. Finally in section 4 we con- 
clude, summarize our main results and propose possible 
extensions of this work. 



II. GEODESICS IN EXPANDING 
SCHWARZSCHILD METRIC 

In order to investigate the effects of expansion on local 
bound systems we must find the geodesies of test particles 
in the appropriate metric that describes the spacetime in 
the vicinity of a point mass M placed in an expanding 
background. Such a metric should interpolate between 
a static Schwarzschild metric at small distances from M 
and a time dependent Friedmann spacetime at large dis- 
tances. The detailed form of this interpolation is not 
unique and there are different a ppr oaches to this prob- 
lem in the literature HH H |2l| • In the Newtonian 
limit (weak field, low velocities) such an interpolating 
metric takes the form: 

ds 2 = (l_^M).^2_ a (i)2. (d(0 2 +p 2. ( ^2 + sm 2^2 )) 



a(t)p- 

where p is the comoving radial coordinate. Using 
r = a(t) ■ p 



(2.1) 



(2.2) 



the geodesies corresponding to the line element (I2.1|) take 
the form 



and 



... a . GM ., n 
-(r r) — +rip=Q 



cp = L 



(2.3) 



(2.4) 



where L is the constant angular momentum per unit 
mass. Therefore the radial equation of motion for a test 
particle in the Newtonian limit considered is 



a 

r = —r 

a r 



GM 



(2.5) 



The same equation of motion is obtained in the Newto- 
nian limit by other interpolations even though the details 
in other limits may vary [Tsj . 
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There is another simple and intuitive (but not rig- 
orous) way to derive the same equation of motion by 
using Gauss's law for gravity with gravitational sources 
the mass M and the integral of the homogeneous source 
p + 3p. This approach leads to 



GM AnG , 
— o~r — {p + 3p)r 



(2.6) 



Using now the Friedmann equation this reduces to 
equation 12.511 . 

Therefore the dynamics of a subluminal test particle 
bound in the gravitational field of a mass M (or equiv- 
alently a two body bound system) in an expanding uni- 
verse can be described by the geodesic equation of motion 
(12.5(1 . In what follows we will study the implications of 
this equation for two body bound systems in various cos- 
mologies. 

As a warm up exercise let us consider the evolution 
of a bound system in an expanding universe with scale 
factor 



a(t) ~ t a 



(2.7) 



where a = const. Let us assume that at some initial time 
to the test particle is at circular orbit with radius tq and 
(p(to) = u>o = ~~?r~ • Then the equation of motion (|2.5() 
may be written in dimcnsionless form as 



a(a — 1) 







(2.8) 



where r = j-, u>q = Wo^o and i = ^. In what follows we 
will omit the bar (777) for convenience but we shall work in 
dimensionless form. Typically for gravitationally bound 
systems in the universe and cosmological timescales we 
have 



GM , 
"T^o » 1 



(2.9) 



(eg for galaxies loq ~ 0(100)) and we may therefore con- 
sider the last term in equation f2.j|j) as a perturbation 
perturbing the circular orbit r$ = 1 to r = 1 + 5r(t). It 
is then easy to show that (see also [Tsj l 



Sr(t) = 



a{a — 1) 
7T572 



(2.10) 



The radius of the orbit tends to increase for a £ (0, 1) 
(decelerating expansion) while the perturbation Sr is neg- 
ative because the homogeneous attractive gravitational 
source (p + 3p > in this decelerating case) within a 
sphere of radius r decreases with time. For a > 1 (accel- 
erating expansion) the perturbation Sr is positive but it 
is decreasing with time because the homogeneous repul- 
sive gravitational source (p + 3>p < in this accelerating 
case) within a sphere of radius r decreases with time. 

To better understand physically the behavior of the 
perturbed orbit we use the Friedmann equation 



This equation shows that the time dependent perturb- 
ing term of equation l|2.8f) comes from the homogeneous 
gravitational source S = ^-(p + 3p)r^ within the unper- 
turbed radius rg. For an equation of state p = wp we 
have 



P 



a -3{l+ W ) „ t -2 



(2.11) 
(2.12) 



The case w > — | co rresponds to a € (0, 1) and the radius 
perturbation (|2.1UI) is negative and decreases with time. 
This is to be expected because the gravitational source 
S ~ p + 3p is positive (attractive) and decreases with 
time (equation l(2.12|l ). Thus Sr is negative (attractive 
perturbation) but decreases with time (the energy den- 
sity decreases with time and so does the homogeneous 
gravitational source). For — l<w<— 4 we have a > 1 
and the radius perturbation I2.10fl is positive decreasing 
with time. This is understood because the gravitational 
source S is negative (repulsive) and decreases with time. 
Thus Sr is positive (repulsive perturbation) but decreases 
with time. Finally for w < —1 we have a < 0. In this 
case the universe has an expanding phase for t < 0, a 
singularity (Big Rip) at t = and a contracting phase 
at t > 0. Our present expanding phase is identified with 
the first phase (t < 0) and the radius perturbation is 
positive and increasing with time. This is to be expected 
because the gravitational source S is negative (repulsive) 
but also increases with time as t approaches the singu- 
larity at t = (see equation I2.12|) . Thus Sr is positive 
(repulsive perturbation) and increases with time. 

This perturbative approach shows that the bound sys- 
tem radial perturbation increases only in a phantom 
(w < — 1, a < 0) expanding background but it can give 
no hint on whether a bound system will eventually be- 
come unbound or not. A non-perturbative approach is 
required to address this question in detail. Such an ap- 
proach should lead to the derivation of the full time evo- 
lution of the radius of a two-body bound system. As 
discussed in the next section this can be achieved by ei- 
ther explicitly solving the equation of motion (|2.5|l or by 
finding the time evolution of the minimum of the effective 
potential derived from equation Q2.8[l . 

Before proceeding with this more appropriate approach 
we will improve on the estimate of the scale factor evo- 
lution by considering a two-component (instead of one- 
component) fluid in a Friedmann universe. Assuming 
the fluid components to be matter with energy density 
p m and dark energy with density p x and equation of state 
p x = wp x we may write the Friedmann equations as [Toj 



8vrG 



[p m+Px ] = H^QKao/af+nKao/a) 3 ^] 

(2.13) 
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and 



AttG 
AttG 



[Pm + Px(l + 3t0)] = 



' Px 



AttG 
l + 3w 



Px 



3w] 



The transition from matter (decelerating) to dark energy 
(accelerating) phase occurs at the transition redshift 



1 -|- z m 



«0 



-(3w + i)n» 



0° 



(2.14) 



At z > z m the solution of the Friedmann equation is 



(2.15) 



where £ is a constant. In the accelerating phase z < z ri 
the corresponding solution is approximated by 



a 



3(l+tu)/2 



(*) 



3 

2 Xt + C 



(2.16) 



where c is an integration constant and x = Ho ^/f2°.. Ex- 
tending the expressions (|2.15() and l|2.16l) up to t TO and 
matching for the evaluation of the integration constants 
£ and c leads to the expression for the scale factor valid 
for t > t m {z < z m ) 



a(t) 



a{t m ) 



-w + (1 + w)t/t m ] ^r+™i 



for t > t m (2.17) 



For phantom energy (w < —1) the scale factor diverges 
at a finite time 



-tm > 



W 



(2.18) 



-3(1+™) 



leading to the Big Rip singularity. Since p x > 
it is easy to see that the phantom energy also diverges at 
t * as 



P*(t) 



p{t m ) 



-W+(l+w)t/t m } 2 



(2.19) 



Using the expression (|2.17|) for the scale factor in the 
two-component fluid universe we are in position to study 
in detail the evolution of bound orbits in phantom and 
quintessence cosmologies. This task will be undertaken 
in the following section. 



As in equation i|2.8[l this equation has been made di- 
mensionless by setting — - 



and j- 

dependent last term proportional to 



t. The time 



a AttG. „ . 2 (l + 3w) 

expresses the gravitating effects of the dark energy in- 
cluded within a sphere of radius r while the dimensionless 
parameter u>q is defined as 



GM 



(3.3) 



We shall use the dimensionless equation of motion i|3.1|) 
with initial conditions r(t = 1) = r m i n and r(t — 1) ~ 
(j"min is the minimum of the effective potential at t = 1 ) 
to study the evolution of the radius of a two body bound 
system which is initially at circular orbit. It will be shown 
that in the phantom case [w < —1) the increasing repul- 
sive effects of the time dependent term of equation <|3.1|) 
lead to a dissociation of the bound system at a critical 
time t r ip. A simplified qualitative approach to this ques- 
tion was made by CKW where no reference was made to 
the equation of motion as it was assumed that the dis- 
sociation occurs when the time dependent dark energy 
gravitational source within the initial radius balances the 

attractive gravitational source ^ . This approach is 
only qualitative and in many cases can lead to incorrect 
results for two reasons: 

• It does not take into account the centrifugal force 

2 

-£q- 

• It implicitly assumes that the radius of the system 
just before dissociation is the same as the initial 
radius r rnin . 

Nevertheless for comparison with our later quantitative 
exact result we will rederive the result of CKW. The bal- 
ance condition of CKW may be written as 



2 27T 2 (l + 3u>) 

w = \TF> - ~n~, 



T' 9 (-w + (1 + u/)t) s 



which leads to 



iy2|l + 3w| 
67HI + w\ 



(3.4) 



(3.5) 



III. BOUND SYSTEM EVOLUTION 

Substituting the derived form of the scale factor in 
equation (|2.17fl to the equation of motion (|2.8|) of a two- 
body gravitating system in an expanding universe we ob- 
tain 



(1 + 3w)r 



9(-w + (l + w)ty 







(3.1) 



where t» is the Big Rip singularity time given by l|2.18|l . 
This is the result of CKW to be compared with our quan- 
titative result derived in what follows. 

The time-dependent effective potential that determines 
the dynamics of the bound system is easily derived from 
equation 13. 1[) to be 



1 



eff 



2r 2 



A(t) 



2„2 



(3.6) 
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where 



A(t) 



y/2|l + 3w| 
3(-iu + (1 + w)t) 



(3.7) 



with w < — 1. At t = 1 the system is assumed to be in 
circular orbit with radius given by the minimum r m ;„ (t) 
of the effective potential of equation Q3.fi [I . The location 
of r mt „ (i) is time dependent and approximates the radius 
of the system at any given time. It is the solution of the 
equation 



l(t) r min — r r, 



- 1 



where 



(3.8) 



(3.9) 



It may be shown (using eg Mathematica [22j |) that this 
equation has a solution only for 



q(tf < = q c 

HK ' ~ 256 y 



(3.10) 



Therefore the time t r i P when the minimum of the poten- 
tial l|3.6[) disappears and the system becomes unbound is 
given by the solution of the equation 



q(t r 



27 
256 



(3.11) 



It is straightforward to solve equation (|3.11|) for t r i P and 
find 



16\/3r v /2|TT3^ 
_ 9 6tt|1 + u;| 



(3.12) 



This result differs from the corresponding result of CKW 
by the factor ~ 3. 

In order to test this result numerically we must solve 
the equation of motion (|3.1I) for a specific bound sys- 
tem, obtain numerically the orbits r(t) and compare the 
analytical prediction for the dissociation time with the 
corresponding time visualized numerically. For concrete- 
ness we consider a phantom cosmology with w = —1.2, 
fl® n = 0.3 and = 0.7. Our goal is to evaluate the 
dimensionless angular velocity lvq corresponding to spe- 
cific systems in the context of this cosmology and use 
it to calculate numerically the corresponding future ra- 
dial evolution. We will consider the scales corresponding 
to three bound gravitational systems: the solar system 
(M = 2 x 10 33 gr, r = 7 x 10 15 cm, u = 3.5 x 10 6 ), the 
Milky Way galaxy (M = 2 x 10 45 gr, r = 5 x 10 22 cra, 
uj = 182) and the Coma Cluster (M = 6 x 10 48 gr, 
ro = 9 x 10 24 cto, ujq — 4.15). In evaluating ojq for the 
above systems we have used the value of t m in the par- 
ticular phantom cosmology considered 



tm 



1.8 x W^h^sec 



(3.13) 
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FIG. 1: The numerical evolution of the radius (continuous 
line) and the analytical evolution of the effective potential 
minimum (dashed line) for the three bound systems. 



We have evaluated the evolution of the radius of the 
above systems using two methods: calculation of the ef- 
fective potential minimum (thick dashed lines of Fig. ^) 
using equation l|3.8f) evolved until the minimum disap- 
pears and explicit numerical evolution of the equation of 
motion (|3.1() (continuous line) evolved up to the dissocia- 
tion time as obtained by CKW (equation l|3.5|l h The nu- 
merical evolution started at t = t m (t = 1 in the dimen- 
sionless form) with initial orbit radius at the minimum of 
the effective potential and negligible radial velocity cho- 
sen such as to minimize radial oscillations. As seen in Fig. 
n there is very good agreement between the numerical 
evolution of the radius (continuous line) and the analyti- 
cal evolution of the effective potential minimum (dashed 
line) for the three bound systems considered. The value 
of the present time to is also indicated on the time axis 
of Fig. 2] In Fig. [21 we show the evolution of the dimen- 
sionless form of the effective potential around the time 
t r i p when the minimum disappears for a bound system 
corresponding to the Milky Way galaxy. As the repulsive 
term destroys the minimum there is a small increase of 
the location of the minimum and then a sudden disap- 
pearance and dissociation of the system. The difference 
between our quantitative prediction for t r i P (indicated by 
the end of the dashed lines in Fig. and the correspond- 
ing qualitative estimate of CKW (indicated by the end 
of the continuous lines) is more prominent for the Coma 
Cluster and less so for the Milky Way. The actual values 
in years for t* — t r i p are shown in Table 1 along with the 
qualitative prediction of CKW. The corresponding dis- 
sociation times t r ip m units of t m are shown in Table 2 
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FIG. 2: The evolution of the dimensionless form of the effec- 
tive potential around the time trip for the Milky Way galaxy. 



FIG. 3: The evolution of the the system of the Milky Way 
galaxy until the effective potential minimum disappears. 



Table 1: The dissociation times differences i* — trip for 
three bound systems in years as predicted by equations 



EDU and ll53 


). The value w = —1.2 was assumed 


System 


t* t r ip (yrs) 


(i* - t r ip)cKW (yrs) 


Solar System 


1.88- 10 4 


6.11 ■ 10 a 


Milky Way 


3.59 • 10 8 


1.17- 10 8 


Coma Cluster 


1.58- 10 1U 


5.14- 10 9 



For the case w = —1.5 considered by CKW, the corre- 
sponding predictions for the Milky Way dissociation are 



rip 



IQQMyrs and (t* — t r i p )cKW — 54Myrs. 



The dissociation times t r i p for the three 



Table 2 

bound systems in units of t m . For w 
t m ~ 5.65h~ 1 Gyrs. 



-1.2, 



System 


trip 1 tm 


Solar System 


6.00 


Milky Way 


5.94 


Coma Cluster 


3.19 



Using the radial equation of motion (|3.1|l along with 
the conservation of angular momentum 



r 2 ip 



(3.14) 



it is straightforward to obtain numerically the full trajec- 
tory corresponding to the evolution of the three bound 
systems and visualize the dissociation process. This dis- 
sociation is demonstrated in Figs. |21and0]for a two body 
bound system corresponding to the Milky Way galaxy. 
In particular, in Fig. |3]the system is evolved until the 
effective potential minimum disappears while the evolu- 
tion in Fig. ^corresponds to the same system but lasts 
until the phantom energy repulsive gravitational force 
balances the attractive gravity of bound matter. Clearly 



the dissociation time is the evolution time of Fig. Awhile 
the evolution in Fig. ^continuous well after the disso- 
ciation as expected based on our analysis. Notice the 
radial (instead of tangential) motion followed after dis- 
sociation which is due to the dominant repulsive gravity 
of phantom energy. 




FIG. 4: The same system as in the previous figure but the 
evolution lasts until the phantom energy balances the attrac- 
tive gravity of bound matter. 

The numerical evolution shown in Figs. ^ an< l 121 nas 
assumed a phantom cosmology with w = —1.2. The dis- 
sociation time t r ip however is sensitive on the value of 
w as shown in equation (|3.12|) . To demonstrate this de- 
pendence Fig. [S] shows the dependence of the relative rip 
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10 - Relative Rip Time Difference 
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+J 

i 

1 - Potential Minimum Disappears 
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FIG. 5: The dependence on w of the relative rip time differ- 
ence. 



time difference defined as 



T 



(3.15) 



(where T is the rotation period of the bound system) on 
w for w < — 1. The continuous curve defines t r i p as the 
time when the effective potential minimum disappears 
while the corresponding definition for the dashed curve 
is the time when the total gravity force vanishes as in 
CKW. Clearly the two curves differ significantly and the 
difference becomes more pronounced as w approaches the 
value w = — 1. 

Our discussion so far was based on the assumption of 
phantom cosmologies (w < —1). This has been due to 
the fact that for w > — 1 (quintessence) the dark en- 
ergy density decreases with time and can therefore not 
destroy the effective potential minimum. Quintessence 
can only cause negligible evolution (decrease) of the ra- 
dius of bound systems. This evolution can be easily ob- 
tained by the perturbative treatment of section 2 by set- 
ting a — 3 j +1 , . To demonstrate the negligible effect of 
expansion on bound systems for non-phantom cosmolo- 
gies we have plotted the effective potential for a Milky 
Way scale system with w — —0.9 at times t = t m and 
t = 3to (Fig. [SJ. The corresponding plot for a Coma 
Cluster scale system is shown in Fig. [7] The decrease of 
the radius is in both cases minor but in the Coma Cluster 
case where ujq is of O(l) it is somewhat more prominent 
as expected from the perturbative result (|2.10[) . 

Finally it is of some interest to investigate the evolu- 
tion of bound systems with planar geometry. Consider a 
test particle at a distance h from a surface with surface 
density a in an expanding universe background. Using 
Gauss's law for simplicity it is straightforward to show 
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FIG. 6: The effective potential for the Milky Way galaxy with 
w = —0.9 at times t — t m and t — 3£o- 
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FIG. 7: The effective potential for the Coma Cluster with 
w = —0.9 at times t = t m and t = 3to- 



that the equation of motion for the particle is 



h = -2ttG(t + --h 
2 a 



(3.16) 



Using equation (|2.17|) for the scale factor in a phantom 
cosmology this may be written as 



t m h = —h + 0h 
where the derivative is with respect to j-, 



h a = 2-KGatt 



and 



3(-w + (l + w)) 2 



(3.17) 



(3.18) 



(3.19) 



8 



Dividing by h and setting 



' h„ 



h, equation hd(t r ip) that dissociates at the time t r i p is found by solv- 



it 17|1 may be written in dimensionless form as 

h=-l+/3h (3.20) 

The effective potential corresponding to this equation has 
obviously no centrifugal term and differs significantly 
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FIG. 8: The effective potential corresponding to bound sys- 
tems with planar geometry for t = t m , 3io and w — —1.2. 
There is reflection symmetry of the potential with respect to 
h = 0. 

from the corresponding effective potential of the spheri- 
cally symmetric case. It is of the form 



r*/~i«»-±>" 



(3.21) 



and is shown schematically in Fig. |SJ It is a reversed har- 
monic oscillator with a time-dependent unstable equilib- 
rium point at h eq = For quintessence (3{t) decreases 
with time and therefore the location of the equilibrium 
point increases with time. Points initially on the left of 
the equilibrium point will remain bound on the attrac- 
tive side of the potential. Points initially on the right 
of the equilibrium (unbound) may also eventually end 
up bound on the left side of the equilibrium point. For 
phantom cosmologies (3(t) increases with time. This im- 
plies that the unstable equilibrium scale h eq (t) decreases 
with time and therefore all scales will eventually become 
larger than h eq (t) and dissociate. The dimensionless scale 



ing the equation 

ha = 



with solution 

U 



3(-U) + t r ip(l + w)) 2 

l + 3w 



T y/2\l + 3w\ 



(3.22) 



(3.23) 



" P ~8V3 \l+w\ 
where T = 4\ / 2hd is the period of the oscillating test 
mass. Notice the similarity of this result with the corre- 
sponding result found for spherically symmetric systems 
(l3~T2l . 



IV. CONCLUSION-OUTLOOK 

We have studied the evolution of bound systems in 
expanding backgrounds and focused on the case of accel- 
erated expansion powered by phantom energy (w < — 1). 
We have found the radial time dependence of bound sys- 
tems in phantom cosmologies and determined the time 
when these systems dissociate due to the repulsive ef- 
fects of phantom energy, as a function of the equation of 
state parameter w. A universal behavior was found for 
the dissociation time for different geometries of bound 
systems. We have also plotted the bound system trajec- 
tory around the time of dissociation and demonstrated 
that the bound systems explode radially outward after 
dissociation. Our results were compared with previous 
corresponding results in the literature and were found to 
be in qualitative but not in quantitative agreement. 

In the present study we have assumed a constant equa- 
tion of state parameter w. The extension of our results 
to the case of a redshift dependent w (w(z)) is straight- 
forward and consists a potentially interesting extension 
of this work. A potential improvement to the accuracy 
of our results may come by using a more accurate metric 
for the interpolation between the Schwarzschild and the 
Friedmann metric [Tsj . Such an improvement would be 
more important for strongly bound high velocity systems. 

The Mathematical^ file used for the production of 
the figures of the paper can be downloaded from |2^] or 
sent by e-mail upon request. 
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